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Abstract 

We examine galactic chemical evolution models which reproduce the present-day and pre- 
solar values of deuterium starting with a primordial value which is consistent with a baryon- 
to-photon ratio of 3 x 10"^*^. We consider various galactic chemical evolution models to 
determine the viability of significant deuterium destruction and which provide a consistent 
age of the galaxy at the time of the formation of the solar system and consequently its 
present day age from nuclear chronometers. These models generally require some amount 
of infall which we take with rates proportional to the gas mass as well as exponentially 
decreasing rates and some initial disk enrichment which we limit to the range of 0% to 30%. 
We present those models which give the observed pre-solar value and present-day value of 
D/H and which lead to a present-day gas fraction of cr = .05 — .2. These models result 
in a braod range for the age of galaxy between 9.8 — 21.6 Gyrs. 



1 Introduction 



One of the successes of the standard big bang model is its abihty to reproduce the observed 
abundances of the hght elements. In two cases, ^He and ^Li, the primordial value can be 
directly determined from observation. The primordial value of ^He is inferred from low 
metaUicity HII regions (see e.g. Pagel et al. 1992; Skillman et al. 1994; Olive & Steigman 
1994). The consistent abundance of "^Li observed in old halo dwarfs is generally assumed 
to be the primordial value (see eg. Spite & Spite, 1993). The primordial abundance of D, 
however, is much harder to test observationally since D is destroyed in stars. A pre-solar 
value determined from meteorites, solar wind studies, and the Jovian atmosphere (see e.g. 
Geiss 1993) and a interstellar value determined from Space Telescope observations (Linsky 
et al. 1992) are all that are known. There is a reported detection of D in a high redshift, 
low metaUicity quasar absorption system (Songaila et al. 1994; Carswell et al. 1994) with 
an abundance which may be the primordial one. Due to the still some what preliminary 
status of this observation and the fact that it can also be interpreted as a if detection in 
which the absorber is displaced in velocity by 80 km s~^ with respect to the quasar (see 
also Steigman 1994; Linsky 1994), we do not fix the primordial abundance with that value 
(we will comment further on this possibility below). The primordial value must therefore 
be explored in the context of chemical evolution models and the observed pre-solar and 
interstellar values (Audouze and Tinsley 1974). 

Several investigations have explored the relationship between chemical evolution and the 
primordial D abundance. (Vidal-Madjar & Gry, 1984; Clayton 1985; Delbourgo-Salvador 
et al., 1985; Vangioni-Flam and Audouze, 1988). It is clear that the degree to which D is 
destroyed will be very sensitive to the chemical evolution model. More recently, Steigman 
and Tosi (1992) explored the acceptable range for the primordial D abundance without 
overproduction of ^He with a given set of chemical evolution models (Tosi, 1988). They 
found that typically deuterium was destroyed by no more than a factor of ~ 2, leading to 
the conclusion that the baryon-to-photon ratio, rj was relatively high 77 >, 4 x 10'^^. A 
different approach was taken by Vangioni-Flam, Olive, & Prantzos (1994), (hereafter VOP). 
They considered a wider range of chemical evolution models starting with a primordial value 
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of D of 7.5 X 10~^ consistent with a baryon-to-photon ratio of 3 x 10~^° a value which is 
more central with respect to '^He and Li consistency. They explored various initial mass 
functions (IMFs) and star formation rates (SFRs) to find chemical histories which could 
reproduce a D destruction factor of ~ 5 to match the presolar value of 2.6 xlO~^ (Geiss 
1993) and the present ISM value of 1.65 xlO~^ (Linsky et al. 1992). They did not, however, 
consider infall in their models and argued against infall since the observations indicate that 
the D abundance has decreased since the pre-solar value. The impact of infall on these types 
of models has not formally been tested. 

A further complication lies in determining the time at which to compare the observed 
abundances with the results from chemical evolution models. This complication is due to the 
large uncertainty in the age of the galaxy. The age of the galaxy and pre-solar epochs can be 
constrained, however, in chemical evolution models which consider, the ratios of long-lived 
radioactive nuclei such as ^^^T/i, ^'^^f/, and ^38^ (Cowan et al. 1987). The advantage to 
including this constraint in with chemical evolution models is two-fold. First, only models 
which can reproduce the ratios observed for these elements are considered viable. Secondly, 
these models give times for the formation of the solar system and today at which direct 
comparisons to observations can be made. 

We will consider the evolution of D using the nuclear chronometer constraint; that is, we 
will only consider models which are capable of producing the correct ratios of ^^^T/i, "^"^^JJ to 
^^^f/ at the same time. Solar system studies indicate a pre-solar D abundance (Geiss 1993) in 
the range of (2.6±1.0) x 10^^. The D abundance has also been determined in the ISM. Lyman 
absorption spectra in nearby stars indicate a present-day D abundance of (.5 — 2.0) x 10~^ 
(see e.g. Vidal-Madjar 1991, Ferlet 1992). A more high precision measurement of D/ H was 
determined recently using the HST (Linsky et al. 1992) which indicates a D abundance of 
1.651^8. For this work, we will adopt the HST measurement for the present-day value. 

Our goal therefore is to explore chemical evolution models which can account for the 
deuterium abundances observed at the pre-solar epoch and today from a primordial value of 
~ 7.5 X 10"^ constrained by nuclear chronometers. It has been argued previously (Reeves 
1991; Reeves 1994) that because of the uranium chronometers, the average galactic nuclear 
activity today is within a factor of two of its average value over the lifetime of the galaxy and 
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therefore may limit the degree to which D is astrated. To test this hypothesis, we consider 
various initial mass functions (IMFs), star formation rates (SFRs) and infall rates which 
in fact can achieve this goal. We also demand that the models reproduce a gas fraction, 
a, in the range of .05 to .2 and require a maximum initial enrichment of 30% with respect 
to present abundances of heavy elements (for a review see Rana, 1991). Finally, we will 
consider the oxygen abundances predicted from such models as a gauge for heavy element 
production. 

2 Chemical Evolution Models 

2.1 General Equations 

In order to examine chemical histories which lead to the observed abundances of D, we 
consider a variety of numerical chemical evolution models. We employ a chemical evolution 
code which avoids the instantaneous recycling approximation (IRA), that is, the actual stellar 
lifetimes are considered in the calculation. We consider a variety of IMFs, SFRs, and infall 
rates. For simplicity, we neglect radial gas flows within the disk for our models, and we do 
not assume that the galactic disk is homogeneous, but restrict our calculation to the solar 
neighborhood. We make no assumptions on the age of the disk. The age is constrained in 
our model by the nuclear chronometers. 

Age determinations, which consider the r-process chronometer pairs 235^^238 ^^^^ 232rpf^^238jj 
have been considered by a number of authors (for review see Cowan et al., 1991). These 
calculations start from the basic equation for gas mass within the galactic disk 

e(t)-^(t) + /(t), (1) 

- mR)i){t - Tm)<j){m) dm. (2) 

In these equations, f{t) is the infall rate onto the disk, 0(m) is the IMF, and il){t) is the 
SFR. rUupp is the upper mass limit on 0(m). is the lifetime of a star of mass m, and rriR 



dMg 
dt ~ 

where e{t) is the ejection rate. 



e{t) = j {m - 



m{t) 
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is the remnant mass of a star of mass m. m(t) is is the mass of a star which at time t is 
returning gas back into the ISM. tur is taken to be (Iben & Tutukov, 1984), 



mfl = .11m + A5Mq m < 6.SMq 

ruR = I.SMq m > G.SM©. (3) 

Stellar lifetimes are determined as in Scalo (1986). 

Equation (1) can be extended to determine the rate of change in the number of nuclear 
species A, 

In this equation, Pa is the number of newly synthesized nuclei of species A per unit mass 
going into star formation. The relative production ratios, P232/-P238, and -P235/-P2.38- arc taken 
to be P232/-P238 = 1-60, and P235/-P238 = 1-16, adopted from Cowan et al. (1987). A^i is the 
rate of decay of nuclear species A. For simplicity, the amount of species "A" which decays 
while locked up in stars has not been considered. This is a good approximation since only the 
largest mass stars will significantly contribute to the amount of "A" produced and returned 
to the ISM and these stars have lifetimes which are short compared with the decay rate of 
^^^Th, '^^^U, and "^^^U. The decay rates used in the calculation are based on the half-lives of 
the respective parent nuclei. For this work, we have adopted half lives of T232 = 1-41 x 10^° 
yrs, r235 = 7.04 x 10^ yrs, and T238 = 4.46 x 10^ yrs. We have also assumed the metaUicity 
of the infall gas, Zf, to be zero. 

Equation (4) can now be solved numerically simultaneously with equation (1). As a 
boundary condition for equation (4), Na{0) is allowed to be non-zero to satisfy observational 
constraints on the metaUicity distribution and age-metallicity relation of stars (Pagel and 
Patchett 1975; Twarog 1980). This is equivalent to assuming an initial burst of star for- 
mation. A^a(O) is a free parameter in our calculation and the value of So, the initial disk 
enrichment, that results from the choice of Na{0) is also calculated. Initial disk enrichments 
of 5*0 ^ .3 are considered acceptable. This corresponds to observational constraints necessary 
to account for the observed metallicities (Lynden-Bell 1975). 
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The final consideration is the abundance of D. We assume no galactic sources of D. 
Furthermore, we consider D to be totally astrated within stars. Under these assumptions, 
the rate of change of D in the disk can be determined by extending equation (1), 

j^iDMG)^-mD + f{t)Df, (5) 

where Df is the mass fraction of D in the infall gas which we assume to be primordial. By 
substitution of equation (1) into this equation, it can be rewritten in the form 

^^-e{t)D + f{t){Df-D). (6) 
This equation is then solved simultaneously with equations (4) and (1). 

2.2 Model Parameters 

The models are divided into three groups based on their SFRs and are closely related to the 
models studied in VOP. In group I, the star formation rate is assumed to be proportional 
to the gas mass (ipit) = vMg)- Group II considers exponential star formation rates = 
z/e"*/"^). We explore decay constants, r, with values r = oo, 10,8,5,3 Gyrs~^. Finally, 
group III models assume a star formation rate proportional to the gas mass squared {ip{t) = 
vMg^)- We allow the constant of proportionality in each of these cases to be a free parameter. 
For most models a simple power-law IMF has been chosen. We do consider the effects of 
varying the IMF. In what we label group lb, we consider an IMF as determined by Tinsley 
(1980). The effect of varying the IMF is also considered for group II models. Group lib 
models employ the IMF given by Scalo (1986) and group Ild models are based on Tinsley 
(1980) with an exponential SFR. The IMF is normalized 

j m(j){m) dm = 1, (7) 

'^low 

where miow and rriupp correspond to the lower and upper mass limits on stars. The values 
of rriiow = -4 and m„pp = 100 have been chosen for this calculation except in the case of 
the Tinsley (1980) and Scalo (1986) IMFs where the value miow = -1 is used. In VOP, a 
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power-law IMF with 0(m) ~ m~^''^ was used and was found to be sufficient for obtaining 
the required D destruction. When we impose the chronometer constraint, we find that in 
many cases not enough D was destroyed and/or the gas mass fraction tended to be too 
low. That is, depending on the assumed SFR, solutions which yield consistent ages, did not 
always destroy enough D and in some cases (eg. models Ila) very few or no solutions were 
found that destroyed D by a total factor of 5 with a sufficiently high gas mass fraction. By 
choosing a flatter IMF, more massive stars which destroy D are produced. Here a power-law 
IMF with (j){m) ~ m~^'^^, has been assumed for most of the models. For model III, an 
even flatter power law ( (j){m) ~ m~^ °^) was required models to destroy enough D to match 
observations. 

Two commonly chosen forms for the infall rate are examined in our calculation. The 
first of these is an infall rate proportional to the gas mass (~ /iMc). Though there is no 
physical reason for assuming this kind of infall rate, under certain assumptions, an infall 
rate proportional to Mq can lead to analytic solutions (see e.g. Clayton 1984,1985). We also 
consider infall rates which arc exponential (~ /xe"*/"^) in models Ic and lie. Decay constants 
of r = oo, 10, 8, and 3 Gyrs^^ are examined. This type of model was previously examined 
by Tosi (1988). A summary of all the models and the various parameters chosen can be 
found in table (1). We have imposed a cut-off on the total disk mass My to be no larger 
than three times the initial disk mass. 

In addition to the constraints outlined above, the abundance of ^^O will be determined 
for each of the model types where the yields of Woosley (1993) have been adopted. For 
models which overproduce oxygen, the value of m„pp, the upper mass limit of stars which 
supernova, is lowered. The impact on the evolution of D by lowering this value to reproduce 
the ^^O abundance will be discussed below. 

2.3 Model Calculation 

The calculation is divided into two parts. In the first, we search for matches for the ratios of 
the nuclear chronometers '^^'^Th/'^^^U and ^^^f/ f^^^U determined from meteoritic abundances. 
The abundance ratios used here are ^^^Th/^^^U = 2.32 and ^^^U/^^^U = .317 (Anders and 
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Ebihara 1982) which have an error of 5 — 10%. In our calculation, we have allowed for an 
error on these values of ±5%. We then look for solutions within this range allowing A^^(O), v 
and /i to vary as free parameters. Chemical histories which produce matches for both ratios 
at the same time are then considered potential solutions. The model time at which the ratios 
match is then considered the pre-solar epoch and subsequently 4.6 Gyrs later, today. Next, 
the D abundance is calculated for models which reproduce these ratios. We test for models 
which can reproduce the D abundance at the times designated as the formation of the solar 
system and "today" and which result in a gas mass fraction a within the range of .05 — .2. 
We examine those models in the next section. 

3 Results 

Most of the model types explored were able to produce chemical evolution scenarios com- 
patible with the observations of the abundance of D. We begin by considering the general 
characteristics that resulted from group I models. Group la models which reproduce the 
observations of D lead to a low of a gas fraction (a ~ .05 — .08) which is only marginally in 
agreement with observations. A higher gas mass fraction can be obtained by further flatten- 
ing the IMF. For example, when 0(m) ~ m^^'^, solutions give a value of o no larger than 
~ 0.05, while for 0(m) ~ m"^'"^ wc find sohition with a as high as 0.2 for model la. It is 
also important to note that our calculational scheme was not designed to maximize the gas 
mass fraction. Instead our sohitions minimize the error in the chronometer abundance ra- 
tios. When this condition is relaxed (though still enforcing ratios within 5% of their observed 
abundances) significantly higher gas mass fractions arc found. 

In group la models, as much as 50% to 70% of the initial disk mass falls onto the disk over 
the galaxy's evolution in these models. In general, the higher the star formation rate, the 
higher the resulting disk mass (see figure la). The solutions shown here have been restricted 
by the requirement that the present D abundance match the ISM measurement to within 
3 %. With a higher star formation rate, more D is destroyed. The only source of D is the 
infall gas. As a result, a larger amount of infall is required to reproduce the observations of 
D for a higher star formation rate. 
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Group lb model solutions closely resemble those of group la. They do result in a slightly 
higher gas fraction than those of la {a <, .09) which is in better agreement with observations. 
The IMF for these models is flatter for the lower mass stars. As a result, less mass is locked 
up in these stars leading to a larger return to the ISM. This allows for more D destruction 
at a higher resulting a. These models require even more infall onto the disk than group la; 
up to 90% of the initial disk mass is required to fall onto the disk to achieve solutions. 

No viable D evolution was found for group Ica models (those with a constant infall rate). 
These models result in a very low o" (cr ~ .03) before enough D is destroyed to match 
observational constraints. Group Icb models can reproduce D observations with very little 
infall (only 27—39% of the initial disk mass), although they do result in a c on the low side 
(o" ~ .05 — .07). Steeper exponents for the infall rate (groups Ice- Ice) can also produce viable 
D evolution scenarios with higher gas fractions at the expense of larger amounts of infall 
required. In order to achieve a gas fraction of o" ~ .1, at least 70 — 90% of the initial disk 
mass is required to fall onto the disk for these model types. Unlike groups la and lb, the 
total mass that results for these model types decreases as a function of the star formation 
rate (as seen in figure la). Note however, that in models such as Ice, there arc solutions 
in which Mt is greater than 3. (They are found at ~ 0.7, where our imposed cut-off is 
visible.) 

Group Ila models which can reproduce the observations of D lead to a gas fractions 
throughout the range of o" ~ .05 — .20. Similar to group I types, the higher the resultant 
SFR, the higher the resultant total mass. Figure lb clearly demonstrates this trend for group 
II models. In addition, as can be seen from figure (2) that in general, as the exponent of the 
star formation rate steepens, we find a lower resulting disk mass and a which themselves are 
correlated within each model type. 

The amount of initial disk enrichment required changes as a function of the decay constant 
of the star formation rate. Group Ilaa solutions (r = oo) result in initial disk enrichments of 
10 — 30%. In general, the steeper the exponent of the SFR, the less initial disk enrichment 
is required. Group Ilae solutions for example all require less than 2% initial disk enrichment. 
Figures 3a and 3b show the range of initial enrichments for solutions as a function of star 
formation rate. Models (except for those with exponential infall rates) with the lowest infall 
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rates and consequently lowest star formation rates lead to the lowest required initial disk 
enrichments. This behavior can be credited to the age at which the pre-solar epoch occurs 
for these solutions. Solutions which have low initial enrichments correspond to later pre- 
solar epochs. Star formation and infall occur for a longer period of time before the pre-solar 
epoch and thus proceed at lower rates but require less of an initial disk enrichment to obtain 
a solution. Figures 4a and 4b show the relationship between the initial disk enrichment, 
5*0, and the pre-solar epoch. The trend of the lower enrichment-later pre-solar epoch is 
demonstrated for the group I and Ila models. Other model groups produce similar results 
(this is true for all models including those with exponential infall rates). From these figures, 
one sees clearly the range of ages for the galaxy at the pre-solar epoch that we obtain in our 
solutions. 

Figures 5a and 5b show a comparison of the evolution of D/H for model groups la, 
lb, Icb and Ilab and their resulting gas fraction. The types are chosen which result in 
a realistic evolution for D/H and satisfy the gas fraction, disk enrichment and nuclear 
chronometers constraints. For comparison, solutions were chosen which had equal pre-solar 
ages. All of these models fall within the errors of both the pre-solar observations and the 
ISM value. There is apparently very little effect on the D/H evolution due to the IMF's 
tested, as these two curves nearly overlap. There is however an effect on the ^^O abundance 
based on the choice of IMF which will be discussed below. These solutions are not unique. 
Other combinations of IMFs, SFRs, infall rates, and resulting pre-solar epochs could produce 
viable D/H evolutions as well. It should be noted that the shape of a solution for the D/H 
evolution for a particular model type does not change for differing pre-solar epochs. The 
general character of each model type remains fixed. 

In all of the cases which produce a realistic D/H evolution, the ^^O abundance is over- 
produced. Adjusting the upper mass limit of stars which supernova {rriupp) for these cases 
can lower the ^^O abundance to its observed levels. Such an adjustment is frequently nec- 
essary when considering exponentially decreasing SFRs (Larson, 1986; Olive, Thieleman, & 
Truran 1987; VOP). The effect of lowering rriupp on the D abundance was also considered 
here. Group I models required lowering the value of rriupp to ~ 20 Mq in order to match the 
solar ^^O abundance. They suffer virtually no change in their predicted D/H evolution by 
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changing this value. In contrast, group Ila models are drastically changed. D is virtually all 
destroyed as a result of lowering rriupp. This result differs from models which assume the in- 
stantaneous recycling approximation (IRA) (see e.g. VOP). Changing rriupp to a lower value 
would result in a lower return fraction. This would mean that less D would be destroyed 
leading to a higher D abundance. In our models, the infall rate has been set proportional 
to the gas mass. This means that less infall will occur in a model where less gas is being 
returned to the ISM due to a lower rriupp- Without this additional source of D, the overall 
effect is to lower the D abundance. Group lib models behaved similarly although the ^^O 
is not overproduced to as great of an extent due to the IMF being normalized more towards 
the low end stars. Group la and lb models have this same assumption. However, in addition, 
the star formation rate is proportional to the gas mass. Thus with less gas returned, the 
star formation rate is lower, leading to less D destruction. This effect compensates for the 
depleted source of D due to a lower infall rate leaving the overall D abundance virtually 
unchanged. Group Ic models, which have exponential infall rates, wind up with an increased 
abundance of D as a result of lowering niupp. In these models, since the infall is independent 
of the gas mass, less stellar processing occurs due to the shift towards the low mass stars 
while D is still supplied by the infall gas. 

Finally, we have considered the effect of changing the D abundance of the infall gas. All 
of the models discussed thus far have considered the D abundance of the infall gas to be 
primordial. In general, it is more difficult to destroy enough D for models which include 
infall. Lowering the abundance of D in the infall gas should, in fact, make it easier to destroy 
D. As a test case, group la models were rerun assuming an infall abundance of D to be 50% 
primordial and also with no D. Figure 6 is a comparison of the D abundance for a group la 
solution. D is destroyed by a factor of ~ 9 for a model with 50% of primordial D abundance 
in the infall gas. With no D in the infall, D is virtually totally destroyed. 

4 Conclusion 

We have shown that D can be destroyed by a factor of ~ 5 under a wide variety of choices 
of SFRs, IMFs, and infall rates even with the additional constraint of nuclear chronometers. 
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It does not seem necessary to constrain the initial SFR to within a factor of 2 of the current 
rate as suggested by Reeves (1991, 1993) in order to reproduce the ratios of the uranium 
chronometers. Though we did find many solutions, we were in fact strongly constrained 
by the chronometers. For example, as in Cowen, Thielemann, & Truran (1987), we could 
not find solutions without infall or without some initial enrichment. The necessity of infall 
was cause for concern with regard to finding solutions with sufficient D destruction factors 
(VOP). However,we did in fact find that significant infall rates can still lead to a decrease 
in the D abundance from the pre-solar value to today. 

In our study, we have always chosen a primordial value oi D / H = 7.5 x 10^^ corresponding 
to a baryon-to photon ratio of ~ 3 x 10~^°. We then selected solutions which yield a present 
deuterium abundance of ~ 1.6 x 10~^ in order to match the ISM measurement (Linsky ct 
al., 1992). In point of fact, we had solutions with both more and less deuterium destruction. 
Thus, our solutions should not be viewed as being tied to our choice of rj or to our assumed 
primordial abundance of D/H. Had we chosen the much higher value oi D/H ^ 2 x 10~^ 
(Songaila et al. 1994; Carswell et al. 1994) we would have still found solutions albeit many 
fewer of them. 

We have also found that adjusting the value of rriupp to match the solar ^^O abundance 
can have a dramatic effect on the evolution of D. In the case of the group Ila models, D 
is destroyed almost entirely by lowering m^pp. It is clear that such solutions exist as they 
require nominally less D destruction than those we have already found. Those models which 
have exponential infall rates (groups Ic, Ild) result in the opposite problem. Not enough D 
is destroyed for these model types when lowering m^pp. 

Any number of our models could be considered good candidates for the evolution of D/H. 
Usually, those models which match observations and result in a gas fraction of o" ~ .1 — .2 
require too much infall onto the disk. Those models with a lower infall rate result in a 
a < .1. (Bear in mind however, our earlier remarks concerning the relationship between the 
IMF and a as well as the effect of our calculational scheme on cr.) These models all assumed 
infall gas with a primordial D abundance. Lowering the abundance of D in the infall gas 
could result in models which can better fit cr without excess infall. 

Our models resulted in the age of the galaxy in the range of 9.8 — 21.6 Gyrs. While 
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nuclear chronometers may not constrain the age of the galaxy very tightly, they have proven 
to be a useful tool in determining viable chemical histories for the reproduction of D obser- 
vations. They give epochs with which to compare with the observations within the models 
and eliminate those combinations of parameters which can not reproduce the observed ratios 
of nuclear chronometers. 
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Table 1: Summary of Model Parameters 
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Figure Captions 



Figure la: Relationship between tiie Total Mass of the galaxy in units of initial disk mass 
and the initial star formation rate for group la (filled circles), Icb (filled boxes), and Ice 
(Filled triangles) model solutions. 

Figure lb: Same as figure (la) for models Ilaa (x's), Ilac (open triangles). Had (open 

boxes), and Ilae (open circles). 
Figure 2: Relationship between the Gas Mass Fraction (a) and the Total Mass of the galaxy 

in units of initial disk mass for groups Ilaa (x's), Ilac (open triangles). Had (open boxes), 

and Ilae (open circles) model solutions. 

Figure 3a: Relationship between the Initial Enrichment (So) and the initial star formation 
rate for models la (filled circles), Icb (filled boxes), and Ice (Filled triangles). The cutoff 
at u — 1 only represents the range of parameters tested. 

Figure 3b: Same as figure (3a) for models Ilaa (x's), Ilac (open triangles). Had (open 

boxes), and Ilae (open circles). 
Figure 4a: Relationship between the Initial Enrichment (5*0) and the resulting age of the 

galaxy at the pre-solar epoch for models la (filled circles), Icb (filled boxes), and Ice 

(Filled triangles). 

Figure 4b: Same as figure (4a) for models Ilaa (x's), Ilac (open triangles). Had (open 
boxes), and Ilae (open circles). 

Figure 5a: Corresponding D evolution for the models in figure (5a) including the observa- 
tions at the pre-solar epoch (Geiss 1993) and today (Linsky et al. 1992). 

Figure 5b: Evolution of the gas fraction for models la (solid line), lb (dotted line), Icb 
(short dashed line), and Ilab (long dashed line) which resulted in a pre-solar epoch of 
9.6Gyrs. 

Figure 6: D evolution for group la model assuming infall D abundances of primordial (solid 
line), 50% primordial (dotted line), and no D (dashed line). 
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